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ABSTRACT 


A  method  for  accelerating:  thin  metallic  plates  to  hypervelocities  has 
been  proposed  by  G.  McCall. In  this  method  a  shock  in  a  propellant 
generates  a  strong  expansion  wave  that  smoothly  accelerates  the  plate. 


We  have  studied  the  hydrodynamics  of  this  process  in  one  dimension,  both 
analytically  and  computationally.  The  metal  was  modeled  as  a  stiffened  gas, 
and  the  corresponding  Riemann  problem  was  solved.  The  asymptotic  behavior  of 
the  solution  was  determined  analytically.  The  one-dimensional  random  choice 
method,  modified  so  that  material  boundaries  are  tracked  and  the  spatial  mesh 
is  refined  locally,  was  used  to  compute  the  flow;  comparison  with  the 
asymptotic  solution  demonstrated  its  accuracy.  With  this  method,  shocks  that 
form  within  the  accelerating  plate  were  accurately  resolved,  so  that  possible 
structural  damage  to  the  plate  could  be  evaluated.  ,  , 
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MODELING  OF  SHOCKLESS  ACCELERATION  OF  THIN  PLATES  USING  A 
TRACKED  RANDOM  CHOICE  METHOD 

* 

Bradley  J.  Plohr 


Introduction 

G.  McCall  has  proposed  a  method  for  accelerating  a  thin  metal  plate  to  exceedingly 
high  velocities  within  a  short  distance.10  These  hypervelocities,  on  the  order  of  1  cm//xs, 
are  difficult  to  achieve  without  destroying  the  plate  since  the  kinetic  energy  of  the  plate 
exceeds  its  vaporization  energy.  The  acceleration  must  be  rapid  enough  to  attain  high 
velocities  within  a  short  distance,  on  the  order  of  10  cm,  but  it  must  be  sufficiently  smooth 
that  strong  shock  waves  do  not  form  inside  the  metal.  Furthermore,  the  compression  of 
the  plate  resulting  from  the  acceleration  must  be  relieved  gradually  so  that  the  plate  does 
not  spall. 

A  schematic  diagram  of  the  experimental  configuration  proposed  by  McCall  is  shown 
in  Fig.  1.  A  tube  containing  a  tungsten  plate  is  partially  filled  with  a  plastic  foam  (CH) 
that  serves  as  a  propellant;  a  low  pressure  cavity  separates  the  foam  from  the  plate.  The 
energy  needed  to  accelerate  the  plate  is  supplied  by  a  planar  shock  passing  through  the 
foam  and  compressing  it.  When  the  shock  reaches  the  end  of  the  foam,  the  foam  expands 
isentropically  into  the  cavity.  At  the  leading  edge  of  this  rarefaction  wav  is  low  density, 
high  velocity  material,  which  is  followed  by  material  with  gradually  increasing  density 
and  decreasing  velocity.  As  succeedingly  denser  portions  of  the  foam  reach  the  plate,  the 
pressure  behind  the  plate  rises,  thereby  gradually  accelerating  the  plate.  Thus  the  low 
pressure  cavity  intervening  between  the  foam  driver  and  the  plate  converts  the  shock  wave 
in  the  foam  into  a  strong  rarefaction  wave  that  smooths  the  acceleration.  Of  course,  the 
rising  pressure  behind  the  plate  produces  compression  waves  in  the  plate  that  in  general 
steepen  into  shock  waves.  McCall  has  estimated  the  formation  time  for  such  shock  waves 
and  has  concluded  that  shockless  acceleration  of  sufficiently  thin  plates  can  be  achieved 
using  his  method. 

A  more  complete  modeling  of  this  hydrodynamic  process  entails  computer  simulation. 
Numerical  methods  for  solving  the  Euler  equations  encounter  three  major  difficulties  in 
modeling  this  flow.  First,  the  mass  density  varies  over  a  vast  range,  from  the  high  density 
in  the  tungsten  plate  and  the  compressed  foam  to  the  near-vacuum  conditions  inside  the 
cavity.  Therefore  Lagrangian  coordinates,  which  place  few  mesh  blocks  in  regions  of  low 
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Figure  1.  Schematic  diagram  of  the  shockless  acceleration  of  a  tungsten  plate  using  a 
plastic  foam  (CH)  propellant  A  planar  shock  in  the  foam  is  incident  from  the  left. 
When  it  reaches  the  end  of  the  foam,  the  foam  expands  into  the  cavity.  The  rarefied 
foam  serves  to  smoothly  accelerate  the  tungsten. 


density,  are  inappropriate  for  the  leading  edge  of  the  rarefaction  wave.  On  the  other 
hand,  the  material  boundary  between  the  metal  and  the  fluid  accelerating  it  is  the  focus 
of  interest.  Whereas  the  position  of  this  boundary  is  naturally  tracked  in  Lagrangian 
coordinates,  it  is  uncertain  to  within  a  mesh  length  on  an  Eulerian  grid;  furthermore, 
mo6t  Eulerian  methods  exhibit  substantial  numerical  diffusion  that  spreads  the  interface 
over  many  mesh  lengths.  Finally,  the  formation  of  shock  waves  within  the  plate  must 
be  accurately  detected,  even  though  the  metal  plate  is  very  thin,  spanning  only  a  few 
computational  mesh  zones. 

In  this  paper  we  have  modeled  the  hydrodynamics  of  McCall’s  problem  using  a  com¬ 
bination  of  methods  that  address  these  difficulties.  The  computation  is  performed  in 
Eulerian  coordinates  to  allow  for  arbitrarily  low  densities  in  the  cavity.  To  correctly  model 
the  material  interfaces,  their  positions  are  tracked  across  the  Eulerian  grid.  Also,  the 
computational  mesh  is  refined  in  the  vicinity  of  the  plate  to  better  resolve  the  flow  in  this 
region.  Finally,  the  underlying  numerical  scheme  for  solving  the  flow  equations,  the  ran¬ 
dom  choice  method,  models  shock  waves  as  sharp  discontinuities  resolved  within  a  mesh 
block,  so  that  shock  formation  can  be  detected  within  the  thin  plate.  This  numerical 
approach  gives  solutions  of  a  much  higher  quality  than  were  previously  available. 

The  plan  of  the  paper  is  as  follows.  In  the  next  section  the  thermodynamic  model  we 
use  for  the  tungsten  plate  is  described.  Then  the  solution  of  the  Riemann  problem,  which  is 
the  key  input  into  the  random  choice  method,  is  described  for  this  thermodynamic  model. 
Next  the  random  choice  method,  as  extended  by  front  tracking  and  mesh  refinement,  is 
explained.  Analytical  results  on  the  short-time  behavior  of  flows  in  which  rarefaction  waves 
impinge  on  a  plate  are  derived  in  the  following  section.  Finally,  the  results  of  numerical 
experiments  for  McCall’s  problem  are  presented. 


Thermodynamics 

The  shockless  acceleration  of  metal  plates  will  be  modeled  as  a  hydrodynamic  process. 
This  requires  the  specification  of  a  suitable  equation  of  state  for  each  component  substance. 
Since  use  of  the  random  choice  method  requires  that  Riemanu  problems  be  solved,  it  is 
advantageous  to  use  an  equation  of  state  that  is  amenable  to  analytic  solution.  The 
propellant  may  be  modeled  as  a  polytropic  gas,  for  which  the  solution  of  the  Riemann 
problem  is  simple,  but  the  metal  plate  cannot  be  adequately  modeled  as  such.  In  lieu 
of  designing  a  Riemann  solver  for  a  general  equation  of  state,  such  as  one  based  on  the 
SESAME  tables,3  we  have  modeled  the  metal  plate  as  a  “stiffened  gas.”  This  equation  of 
state  gives  the  correct  semi-quantitative  behavior  of  the  plate.  In  conjunction  with  a  spall 
model  to  give  yield  strains  and  an  estimate  of  melting  and  boiling  points,  it  can  be  used 
to  asses  the  feasibility  of  shockless  acceleration. 

The  equation  of  state  for  a  stiffened  gas8  is  given  by  the  relation 

E_  1  P  +  7Poo 
7  —  1  P 

between  the  specific  (internal)  energy  E ,  the  pressure  p,  and  the  mass  density  p.  Here  7 
and  poo  are  prescribed  constants  that  characterize  the  gas,  7  being  dimensionless  and  Poo 
being  a  pressure.  In  terms  of  the  Gruneisen  coefficient 


f\9E)t 

we  have  7  =  T  +  1.  (Note,  however,  that  7  is  not  the  adiabatic  exponent  7 s,  which  is 
defined  by  c3  =  7 sp/p\  see  (2.2)  below.) 

The  isentropes  for  a  stiffened  gas  may  be  determined  by  noting  that 


dE  +  p  d 


(;)  = 


P+Poo  1 


p  7-1 

Thus  the  entropy  S  is  most  generally  given  by 


d  log 


(^=). 


for  some  function  E,  with  the  corresponding  temperature  T  being 


(2.1) 


T  = 


P  +  /-(log(^)) 


The  speed  of  sound  c,  defined  by 


P  +  Poo 


Fast  and  slow  rarefaction  waves,  which  correspond  to  the  fast  and  slow  characteristic 
speeds  v  ±  c,  satisfy  the  equations 

dt;  =F  —  =  0 


dE  +  p  d 


By  (2.1)  and  (2.2)  we  determine  that 


dp  _  2 

pc  ”7  —  1 ' 


Therefore 


v  T - c  =  const. 

-7  —  1 


(p  +  Poo)/p'7  =  const. 

across  rarefactions  corresponding  to  the  characteristic  speeds 

v  ±  c  . 

The  Hugoniot  locus  is  given  by  solving  the  equation 


E  -  Eo  +  ~(p  +  po)  ( - ~>)  =  0 

*  \  P  Po J 


for  p  in  terms  of  po,  po,  and  p.  Simple  algebra  yields  the  result 

P  _  (7  + l)(P  +  Poo)  +  (l  -  l)(PO  +  Poo) 


po  n  +  !)(po  +  Poo)  +  h  -  i)(p  +  p<xj) 

As  for  any  equation  of  state,  the  speed  a  of  the  corresponding  fast  and  slow  shocks  may 
be  determined  from  the  relations 


p(v  -  o)  =  Tm  =  po{v0  -  a)  , 


where 


■(- 


P  ~  Po 

i Ip  -  1/po. 


-5- 


is  the  mass  flux  across  the  wave. 

To  demonstrate  the  adequacy  of  using  a  stiffened  gas  equation  of  state  to  model  a 
metal  plate  made  from  tungsten  we  compare  to  experimental  data.1  These  data  show  that 
a  plot  of  shock  speed  U,  =  \a  —  t>o|  vs.  particle  speed  Up  =  \v  —  vol  along  the  principal 
Hugoniot  of  tungsten  is  well  approximated  by  a  straight  line 


U.  w  c0  +  <*UP  (2.9) 

for  pressures  up  to  3  Mbar.  Here  Cq  =  0.404  cm/ps  is  (approximately)  the  sound  speed  in 
tungsten  at  standard  conditions  and  a  =  1.23  is  an  empirical  dimensionless  constant.  For 
a  general  equation  of  state  one  finds  from  (2.7)  and  (2.8)  that 

Up  _  p  -  Po 
U.  ~  p 


U.UP  = 


P~  Po 


Using  the  Hugoniot  relation  (2.6)  for  a  stiffened  gas,  one  can  eliminate  p  and  p  to  find 

U?  =  ^U.u,  +  4 


+  !±ip 

+  4  . 


(2.10) 


For  our  calculations  modeling  tungsten  we  have  used  7  =  3.14  and  =  1.0  Mbar.  A 
plot  comparing  the  empirical  linear  fit  to  the  stiffened  gas  relation  (2.10)  for  pressures  up 
to  3  Mbar  is  shown  in  Fig.  2.  The  stiffened  gas  relation  coincides  with  the  linear  fit  to 
within  3%. 


particle  speed  Up 


Figure  2.  Comparison  of  the  stiffened  gas  equation  of  state  to  experimental  data  for 
tungsten.  The  shock  speed  Us  along  the  principle  Hugoniot  is  plotted  with  respect  to  the 
particle  speed  Up  behind  the  shock.  Speeds  are  measured  in  cm/|is ,  and  the  range  of 
shock  strengths  corresponds  to  overpressures  up  to  3  Mbar.  The  stiffened  gas  is  defined 
by  the  parameter  choices  y=3.14  and  p~=  1.0  Mbar.  The  comparison  shows  a  max¬ 
imum  difference  of  3%. 


Riemann  Problems 


A  Riemann  problem  for  a  system  of  conservation  laws  in  one  space  dimension  is  an 
initial-value  problem  with  scale  invariant  initial  data,  i.e.,  data  at  t  =  to  that  is  constant 
to  either  side  of  a  jump  discontinuity  at  x  =  xo-  We  use  the  subscript  l  to  label  the  initial 
state  on  the  left  of  io,  and  r  for  the  initial  state  on  the  right  of  xo-  The  solution  of  a 
Riemann  problem  is  constant  along  each  ray  £  =  (x—xo)/(t  —  t0)  =  const,  through  (x0,t0); 
thus  we  may  specify  the  solution  as  a  function  of  the  speed  £. 

The  solution  of  a  Riemann  problem  for  a  stiffened  gas  requires  only  a  simple  modi¬ 
fication  of  the  well-known  solution  for  a  polytropic  gas.  It  consists  of  a  slow  (i.e.,  v  —  c) 
rarefaction  or  shock  wave  on  the  left,  a  fast  (i.e.,  t>  +  c)  rarefaction  or  shock  wave  on 
the  right,  and  a  contact  discontinuity  (moving  with  velocity  v)  in  between,  as  indicated 
schematically  in  Fig.  3.  It  is  convenient  to  speak  of  the  left  side  of  a  slow  wave  and  the 
right  side  of  a  fast  wave  as  being  “ahead”  of  the  wave,  the  opposite  sides  being  “behind” 
the  wave;  the  corresponding  states  are  labeled  by  the  subscripts  a  and  6.  A  slow  or  fast 
wave  is  determined  by  the  state  ahead  of  it  together  with  one  other  parameter,  such  as  the 
pressure  behind  it.  For  the  left  (resp.,  right)  wave  in  the  solution  of  a  Riemann  problem, 
the  state  ahead  of  the  wave  is  the  left  (resp.,  right)  initial  state.  Since  the  particle  velocity 
and  the  pressure  is  continuous  across  a  contact  discontinuity,  the  velocities  and  pressures 
behind  the  left  and  right  waves  must  coincide.  Thus  solving  a  Riemann  problem  amounts 
to  finding  the  pressure  in  the  middle  between  the  left  and  right  waves  such  that  the  particle 
velocities  behind  the  waves  are  equal. 

For  a  shock  wave,  conservation  of  momentum  may  be  written 

^m(va  -  vt)  +  pa  ~  Pfc  =  0  •  (3-1) 

Here  m  is  defined  as  a  function  of  the  pressure  pb  behind  the  Bhock  and  the  state  ahead  of 
the  shock  through  equations  (2.6)  and  (2.8),  so  long  as  p&  >  p„.  For  a  rarefaction  wave, 
we  may  define  m  through  the  equation 

2 

-m - -(ca  -  Cfc)  +  Pa  -  Pb  =  0  , 

7~1 

so  that  by  (2.3)  we  find  that  (3.1)  again  holds.  With  this  definition  and  equations  (2.4) 
and  (2.2),  m  may  be  expressed  as  a  function  of  p&  and  the  state  ahead  of  the  wave,  so  long 
as  Pb  <  Pa-  More  explicitly  we  have 


where 


m  = 


{laPa(Pa  4"  Poo, a 


Pb  +  Poo,  b 
Pa  +  Poo, a 


^(r>  7) 


! 


7-1  1  -  r 

2  7  1  -rh-O/^T 

7+2  izlV/2 

27  27  / 


if  r  <  1 
if  r  >  1  . 


(3.2) 
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Figure  3.  The  computational  grid  for  the  random  choice  method.  The  reconstructed 
solution  is  constant  in  the  cells  centered  on  the  circles  and  bounded  by  the  vertical  hash 
marks.  The  solution  of  the  initial  value  problem  with  piecewise  constant  initial  data 
involves  solving  a  sequence  of  Riemann  problems  originating  at  the  boundaries  between 
cells.  The  solution  of  a  Riemann  problem  is  constant  along  rays  through  its  origin  and, 
for  gas  dynamics,  consists  of  a  shock  or  rarefaction  wave  (indicated  by  an  individual  ray 
or  a  fan  of  rays,  respectively)  on  each  side  of  a  contact  discontinuity  (indicated  by  a 
dashed  ray).  The  numerical  solution  after  one  time  step  is  obtained  by  sampling  the 
Riemann  solutions  at  the  points  indicated  by  the  stars,  whose  shift  relative  to  the  cell 
centers  is  random  in  time. 
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Let  v,  and  p,  denote  the  common  values  of  the  particle  velocity  and  pressure  on  the 
two  sides  of  the  contact  discontinuity.  Then  by  (3.1)  we  have 

-mT(vr  -  v,)  +  pr  -  p.  =  0 


+m,(vj  -  v*)  +  pi  -  p,  =  0 

for  the  right  and  left  waves,  where  mr  is  expressed  through  (3.2)  in  terms  of  pb  =  p, 
and  the  density  and  pressure  pr  and  pr  ahead  of  the  wave,  and  similarly  for  m,.  These 
equations  may  be  written 

»l-»,  mIPr+mrPl 

1/mi  +  l/mr  m\  +  mr  '  '  ' 


P/  ~  Pr  +  +  mTvr 

* - ^  TZZ - •  (3-4) 

Therefore  to  solve  a  Riemann  problem  we  first  use  the  Godunov-Chorin6’ 3  fixed-point 
iteration  to  solve  (3.3)  for  the  pressure  p*  between  the  left  and  right  waves.  Then  (3.4) 
gives  the  velocity  v*  in  the  middle,  and  for  both  the  fast  and  the  slow  wave  the  state  along 
a  ray  moving  at  a  specified  speed  {  may  be  determined  using  either  equations  (2.2-5) 
or  equations  (2.6-7),  according  as  the  wave  is  a  rarefaction  or  a  shock.  Except  for  the 
modifications  to  the  formulae  (2.2),  (2.4),  (2.6),  and  (3.2),  the  construction  of  the  solution 
of  the  Riemann  problem  is  identical  to  that  for  a  polytropic  gas  (for  which  p ^  =  0) . 


Tracked  Random  Choice  Method 

The  random  choice  method  is  a  numerical  method  for  solving  the  general  initial-value 
problem  for  a  system  of  conservation  laws.  It  is  based  on  the  constructive  existence  proof 
for  solutions  of  conservation  laws  of  Glimm,4  and  was  adapted  for  numerical  computation 
by  Chorin.3 

In  the  random  choice  method  for  the  hyperbolic  system  of  conservation  laws 

+  /(«)*  =  0  (4.1) 

the  solution  is  determined  at  the  lattice  of  space-time  points  (xy,  tn); 

Xj  =  (xy_1/2  +  xy+i/a)/2 

is  the  center  of  the  mesh  cell  from  xy_,/2  to  xy+i/2.  Given  the  solution  values  u"  = 
u(ij,tn)  at  time  level  »,  the  solution  values  u”+1  at  level  ft  +  1  are  determined  in  three 
stages. 

First,  an  approximation  u"(x)  for  the  solution  at  time  tn  and  an  arbitrary  position  x 
is  reconstructed  by  defining  un  to  be  piecewise  constant  on  mesh  cells: 

un(x)  =  u *  if  Xy-i/2  <  x  <  Zj+i/2  •  (4-2) 

Next,  consider  solving  the  initial-value  problem  for  (4.1)  with  initial  data  (4.2)  at 
time  tn.  For  t  —  tn  sufficiently  small,  the  solution  is  given  in  terms  of  the  solutions  of  a 
sequence  of  Riemann  problems  between  the  states  u"  and  u"+1  and  centered  at  the  points 
(xy+i/2,*n).  The  situation  is  depicted  in  Fig.  3.  In  fact,  the  waves  in  the  solutions  of  the 
Riemann  problems  centered  at  neighboring  points  (xy_i /2,tn)  and  (xy+i/2,fn)  will  not 
interact  so  long  as  the  Courant-Friedrichs-Levy  condition 

t-tn  <  tCFL  =  \  min  (xy+1/2  -  xy_i/2)  / Ay  (4.3) 

i  J 

holds,  where  Ay  is  the  maximum  wave  speed,  in  absolute  value,  in  the  solution  of  the 
Riemann  problem  centered  at  (xy+1/a,tn).  We  have  taken  tn+i  —  tn  =  C  •  <cfl  with  the 
Courant  number  C  kept  between  0.7  and  0.9. 

Finally,  an  equidistributed  random  number  r?n+i  is  chosen  from  the  interval 

-1/2  <  fjn+l  <  1/2  , 

and  the  solution  value  tt"+l  is  taken  to  be  the  value  u*  of  the  solution  of  the  above 
initial-value  problem  sampled  at  the  position 

X.  =  Xj  +  T)n+1  •  (xy+i/2  ~  Xy_j/2) 

at  time  tn.  In  our  computations  we  have  used  the  Van  der  Corput  sequence3  to  obtain  rj. 

The  solution  given  by  the  random  choice  method  has  the  advantage,  when  compared 
to  standard  finite  difference  methods,  of  preserving  the  structure  of  discontinuous  waves  in 


he  .  ulu'ion:  sho^k  waves  and  cin.urt  discontinuities  remain  perfectly  sharp.  The  method 
.  .horeioie  able  to  detect  formation  of  shock  waves  within  regions  extending  only  a  few 
mesh  cells.  This  property  is  crucial  for  our  application  to  shockless  acceleration  of  thin 
plates. 

The  price  paid  for  perfect  resolution  of  discontinuous  waves  is  that  their  positions  are 
not  exactly  correct.  On  the  average  the  position  of  a  discontinuous  wave  in  the  numerical 
solution  is  correct,  but  the  wave  performs  a  random  walk  about  this  correct  position.  To 
remedy  this  in  situations  where  maintaining  the  position  of  certain  prominent  waves  is 
critical,  Glimm,  Marchesin,  and  McBryan  introduced  the  tracked  random  choice  method.6 

In  a  front  tracking  method,  certain  waves  in  the  flow,  typically  strong  discontinuous 
waves,  are  resolved  on  a  subgrid  level.  With  each  tracked  front  is  associated  a  position 
and  a  wave  type.  The  wave  type  determines  how  the  tracked  front  moves  and  interacts 
with  the  fluid.  For  example,  in  gas  dynamics  the  possible  wave  types  include  slow  and  fast 
pressure  waves  and  contact  discontinuities. 

Let  x"  denote  the  position  of  a  particular  tracked  front  at  time  tn.  If  the  tracked 
front  lies  in  the  jth  mesh  cell  at  time  tn,  i.e.,  Xj-1/2  <  x"  <  xJ  +  1/2,  then  the  solution 
values  ,  u” +1,  and  are  obtained  in  a  fashion  different  from  that  just  described. 

The  piecewise  constant  data  un  is  defined  by 


if  Xj-3/2  <x<x? 
if  X?  <  X  <  x;+3/2 

for  x  in  the  interval  xy_3/2  <  x  <  Xy+3/2,  and  by  (4.2)  for  x  outside  this  interval.  As 
before,  the  initial-value  problem  with  initial  data  un  at  time  tn  is  solved  up  to  time  tn+  j. 
The  position  of  the  tracked  front  at  time  tn+l  is 


x)  = 

K- 

luy+i 

<  X 

<  Xj +3/2 

x"+1  =  Xj  +  s  •  (t„  +  i  -  tn)  , 


where  the  velocity  s  is  determined  according  to  its  wave  type.  Finally,  the  solution  value 
Uy^, 1  is  taken  to  be  the  value  it.  of  the  solution  of  this  initial-value  problem  at  time  f„+i 
sampled  at  the  position 

=  xj- 3/2  +  (*7n+l  +  1/2)  (x"+1  -  Xj _  3 / 2 )  . 

Similarly  u"  ^j1  is  obtained  by  sampling  at 

x.  =  x("  +  1  +  (r/n+l  +  1/2)  (xj  +  3/2  x"+1)  • 


The  choice  for  the  value  of  Uy  +  1  seems  somewhat  arbitrary.  For  instance,  one  could 
take  Uy  +  1  to  be  Uy  or  Uy  +  j1  according  as  xy  <  or  Ifl+1  <  x;  ■  have  found 

that  with  this  choice,  however,  a  strong  wave  crossing  the  tracked  wave  does  not  separate 
correctly  from  the  tracked  wave  after  the  interaction.  This  effect  is  rather  dramatic  in 
calculations  where  a  tracked  shock  wave  is  followed  closely  by  a  strong  rarefaction  wave. 
To  remedy  this  we  instead  define  u"+1  to  be  u"^1  or  u"++  1  according  as  x;  <  x”  +  1  or 
x”  T  1  <  x.j,  where  1  anrl  u”4f  1  are  the  states  on  either  side  of  the  tracked  front. 
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In  gas  dynamics  the  possible  wave  types  are:  contact  discontinuities,  moving  with  the 
particle  velocity  v;  slow  and  fast  pressure  waves,  which  either  follow  sound  waves  moving 
with  velocity  v±c  in  case  the  solution  is  continuous  across  the  wave  (the  rarefaction  case), 
or  move  with  the  shock  velocity  a  if  the  wave  is  discontinuous  (the  shock  case);  Neumann 
and  Dirichlet  boundaries,  which  supply  boundary  conditions;  and  movable  walls,  which 
are  ascribed  a  mass  density  and  move  according  to  Newton’s  law.  The  velocity  of  a 
slow  or  fast  pressure  wave  or  contact  discontinuity,  is  calculated  by  solving  the  Riemann 
problem  between  and  centered  at  (x^,fn),  and  determining  the  velocity  of 

the  corresponding  wave  in  this  solution.  (We  track  the  leading  edge  of  a  rarefaction 
wave.)  Neumann  and  Dirichlet  boundaries  are  stationary,  and  their  interaction  with  the 
flow  is  determined  by  solving  a  Riemann  problem  between  the  fluid  state  neighboring 
the  boundary  and,  respectively,  the  same  state  with  a  reflected  velocity  or  the  imposed 
boundary  state.  The  evolution  of  a  movable  wall  is  determined  as  follows. 

Let  the  wall  have  mass  per  unit  area  M  and  velocity  v£a,|  at  time  tn.  Suppose  for 
concreteness  that  the  fluid  lies  on  the  right  side  of  the  wall,  so  that  the  solution  value 
inside  the  fluid  is  ur  =  u"+1.  Define  the  virtual  state  to  have  the  same  pressure  and 
density  as  ur  and  to  have  the  velocity 

vl  =  2tCal,  -  t>r 

obtained  by  reflection  through  the  moving  wall.  Then  solve  the  Riemann  problem  between 
u;  and  ur  centered  at  (x",*,,)  and  sample  it  at  the  position 

X«  =  If  +  vwall  ‘  (^n+i  ~  fit) 

to  obtain  the  state  u,  on  the  fluid  side  of  the  wall  at  time  fn+j.  Because  vj  is  obtained 
by  reflection,  the  velocity  v,  of  the  fluid  at  the  wall  coincides  with  v”all,  so  that  the 
correct  boundary  conditions  have  been  maintained.  The  pressure  in  the  fluid  accelerates 
the  movable  wall  according  to  Newton’s  law.  If  At  =  tn+i  -  tn  is  small  enough  (more 
precisely,  as  may  be  deduced  from  the  analytic  solution,9  if  p.At/Afc.  is  small)  we  may 
approximate  this  acceleration  by 

„«+!  =  „»  .  El±Ei  . 

Vwall  vwall  ^  2  ' 

The  speed  ,s  of  the  moving  wall  during  this  time  interval  is  then  taken  to  be 

*  =  \  Wtil  +  <„.)  • 

To  obtain  acceptable  resolution  of  the  waves  in  the  thin  plate  the  computational  grid 
must  be  refined  in  the  vicinity  of  the  plate.  The  positions  of  the  the  edges  of  the  mesh 
cells  x;  +  1/2,  j  =  0, ...,N  are  defined  as  follows.  The  grid  spans  the  interval  [xm;n,xmax] 
and  has  a  subinterval  [a,  6]  over  which  the  mesh  is  refined  by  the  refinement  ratio  R.  To 
smooth  the  transition  between  the  coarse  and  fine  meshes,  a  molification  parameter  £  is 
also  specified,  with  £  =  0  corresponding  to  a  sharp  transition.  We  define 


_  / 1/* 


■<«  =  {! 


if  0  <  £  <  1 
otherwise 


and  let  pt  be  the  convolution  pt  =  St  *  po  of  Po  with  an  approximate  6-function.  We  have 
used 


5,(f)  =  \  h  (|) 


where 


MO  = 


With  these  definitions,  the  zone  edges  are  given  by 


C/+1/3  —  *mln 


+  (lout  X»,»)  •  £  dt  pt( 0  j  J*"  di  pM) 


where 


&«&+£•  tin  ~  &)  • 


It  remains  to  define  Zo  and  Zn  •  The  transitions  occurring  between  the  coarse  and  fine 
meshes  at  Z  —  0  and  £  =  1  in  (4.4)  are  to  correspond  to  x  =  a  and  x  =  b  through  (4.5). 
For  c  ^  0  this  gives  a  pair  of  nonlinear  equations  to  be  solved  for  Zo  and  Zn  in  terms  of  a 
and  6.  When  6  =  0,  however,  these  equations  reduce  to  linear  equations,  with  the  solution 

1  4  -  Xmin 
~  #  b-  a 


1  xmmx  -  b 

for  simplicity  we  use  these  equations  to  define  Zo  and  Zn ■ 


Analytical  Results 

As  a  preliminary  to  solving  the  shockless  acceleration  problem,  we  study  a  similar 
problem  for  which  analytical  information  is  available.  This  is  the  problem  of  the  reflection 
from  a  rigid  wall  of  a  rarefaction  wave  expanding  into  a  vacuum.  The  short-time  asymptotic 
solution  of  this  problem  has  been  described  by  Greenspan  and  Butler;7  they  reduce  the 
problem  to  the  solution  of  a  boundary-value  problem  for  a  system  of  ordinary  differential 
equations,  which  they  solve  numerically.  These  equations  may,  in  fact,  be  integrated 
analytically,  and  we  describe  this  solution  here. 

The  flow  pattern  is  shown  in  Fig.  4.  A  rarefaction  wave  in  a  polytropic  gas  with 
adiabatic  exponent  -7  connects  the  state  p  —  po ,  p  =  po,  v  =  Vo  to  the  vacuum  p  —  0, 
p  =  0;  it  originates  at  x  =  0  at  time  t  =  0.  The  rarefaction  wave  corresponds  to  the  slow 
(j.e.,  t;  -  c)  family  of  characteristics,  so  within  the  wave  the  fluid  state  is 

p  =  Po(c/co)3/('r_1)  , 

P  =  Pc(</eo)2V('’-,)  •  (5-1) 


v  =  v0  +  - -(Co  -  e)  , 

If  -  1 


where 


v  -  n  • 


by  virtue  of  Kqs.  2.2-  5.  Since  the  tail  of  the  wave  corresponds  to  c 
arrives  at  the  wall  located  at  1  =  i0  at  time 


=  0,  the  rarefaction 


Xo  I  v0  + 


vr 


At  this  time  a  reflected  shock  wave  forms  to  decelerate  the  expanding  fluid. 

Following  Greenspan  and  Butler  we  seek  a  scaling  solution  for  the  flow  between  the 
reflected  shock  and  the  wall  in  terms  of  the  variables  r  -  t  —  t q  and  (  =  (x  -  z0)/(t  to). 
Ahead  of  the  reflected  wave  the  state  is  given  by  Eqs.  5.1,  where  c  may  be  written 

with  £0  -  Zo't o-  Thus  the  state  just  ahead  of  the  reflected  shock  wave  satisfies  pa 
()  ((r/f0)2/l'7  1 1 ) ,  pa  O  ((r  In)3 i)),  and  va  -  0(1)  as  r  *0.  The  state  p  pk. 
p  pk,  v  ~  v,  just  behind  the  reflected  shock  wave  must  satisfy  the  Hugoniot  relation 
Fiq.  2  6  together  with 

(n  Ua)3  :  (pb  pa )  ( 1  /  Pa  1/Pb)  .  (•'»■- 

which  is  a  consequence  of  momentum  conservation.  Assuming  a  non-negligible  part  of  the 
deceleration  of  the  fluid  occurs  across  the  shock  wave,  so  that  its  strength  does  not  vanu-h 
0,  the  left-hand  side  of  this  last  equation  is  0(1).  It  follows  that 


O  ( (r  't0 


2  [i  1) 


as  t  ♦ 


reflected  shock 


Figure  4.  Reflection  of  a  rarefaction  from  a  rigid  wall  A  rarefaction  wave  in  a  polytro 
pic  gas  connects  the  state  with  p  =  po.  p  -  p o,  and  v  =  vn  to  a  vacuum  state  The  rarefac 
tu >n  originates  at  x  = 0 ,  and  its  tail  strikes  a  rigid  wall  at  »  =  to  at  time  t  =  rn,  whereupon 
a  reflected  shock  is  formed. 


and  hence  that 


Pb  =  - - ~Pa  •  (l  +  0(r(t o)2) 

If  -  1 

as  t  — »  0.  Thus  the  reflected  shock  is  infinite  in  strength  at  the  time  of  its  formation. 

In  this  way,  Greenspan  and  Butler  are  led  to  expect  that  the  solution  between  the 
reflected  shock  wave  and  the  wall  is,  asymptotically  as  r  — *  0,  of  the  form 

p=(’/>oy,h-,,M) .  (s.3) 

V  =  0(0  ; 

here  (  lies  in  the  range  a0  <  £  <  0,  Oo  being  the  limiting  velocity  of  the  reflected  shock 
wave  as  r  — ►  0.  By  assuming  that  the  solution  is  of  this  form,  they  find  a  system  of 
nonlinear  ordinary  differential  equations  to  be  satisfied  by  p,  p ,  and  0,  and  they  solve  the 
equations  numerically.  However,  one  can  find  an  analytic  solution  of  these  equations  in 
terms  of  incomplete  beta  functions,  as  we  now  show. 

The  flow  is  governed  by  the  characteristic  equations 

vt  ±  — pi  +  (t>  ±  c)  ( vz  ±  —  pt  ]  --  0  , 
pc  \  pc  J 

Pt  -  c2pt  +  V  {pz  -  C7p1)  0  . 

It  is  subject  to  no-flow  boundary  conditions  t>  =  0  at  £  =  0,  and  to  the  boundary  conditions 


at  ‘  ex.  ■  these  conditions  follow  from  mass  conservation,  Kq  2  7,  and  moment  inn  <<</, 
servation,  Kq.  5.2,  in  the  strong-shock  limit. 

I'pon  substituting  the  polytropic  gas  relation  r2  ip  p ,  the  <  har <u' tenst .<  equation" 
may  be  written 


When  the  solution  is  assumed  to  be  of  the  form  of  Eq.  5.3,  v  =  v  and  e  —  e  -  {ip/p)1^2 
are  independent  of  r,  so  the  last  equation  becomes 


2  ^  1  _  2  c 

i  -  ic  1(1  -  i) 6  -  £ 


therefore  the  previous  pair  of  equations  reads 


(v±i-t)[v,± 


2  c2 

nr(nr  -  l)  »  -  €  ’ 


where  the  prime  denotes  differentiation  with  respect  to  £. 

In  terms  of  e  —  c  and  the  Mach  number  M  =  (v  ~  ()/c,  Eqs.  5.5  are 


M(M  ±1)  1  +  cM‘  +  ±  e'j  = 


l(l  ~  1)  ’ 


cM'--  —  (l  -  2~iiM2)  (l  -  M7)  1 

2i  t 

and 

Me'  =  i(l  -  M2)  ‘  , 

1 

while  Eq.  5.4,  when  combined  with  Eq.  5.7,  becomes 

P 

We  have  introduced  the  parameter  *  —  \(l  ~  1 ) / (n^2  +  1)  for  later  convenience. 
The  ratio  of  Eqs.  5.6  and  5.7  dictates  that 

=  1(1  2-r.M1)  M1, 

a  log  c  « 


which  integrates  to  yield 


c  c* 


\  M2 


2  1  2i(M2  \  * 

2  1  2~t<M?  ) 


Here  ,\ff  (rfc  On )/ch  v  (-7  1  )/2"7  is  the  Mach  number  just  behind  the  reflected 

"ho<  k  wave.  Similarly, 

dM 2  1  1  i  ,  .  w2\ 


d  log  p  2*m 


(1  M2) 


.  .  /  1  2 \ 

;n.  21(m2) 


2  ,  !  ”  «> 


^•VV.V.WWjV^  V’/  b  '.".>7/.  c.v.vy.  '  ■  V  /  .  • 


By  using  Eq.  5.8,  Eq.  5.6  may  be  written 


(1  -  M2)  (1  -  2'1eM2y{l-t){\f2)-{*+t)dXf 2  = 

Consequently,  if  we  define 

V'-r(y)  =  /  (1  -  *)(1  -  2-yc  ■  z)~^1~^z~^+e^dz 
Jo 


(where,  again,  <  =  5(7  -  l)/('72  +  1)),  then 


(5.10) 


Since  is  a  monotone  function,  Eq.  5.10  may  be  solved  for  M2  in  terms  of  £,  which 
gives  c  and  p  as  functions  of  £  through  Eqs.  5.8  and  5.9.  Furthermore,  by  imposing  the 
boundary  condition  that  M  =  0  at  £  =  0,  we  obtain  an  equation  for  00/ £0  as  a  function 

of  -7: 


U»  _  =  _oo_ 

t>o  ■+  2c0/(t  -  l)  £0 


1  7  +  1  1  /l-27cMb2V  i 

v/2nr(nr  -  0  i)  V  W?  ) 


(5.11) 


Notice  that  the  function  tl>~,  may  be  expressed  in  terms  of  incomplete  beta  functions, 
which  are  defined  by 

9;  y)  =  /  «p“ 1  (1  ~  2),_1ds 

Jo 

when  p  and  q  have  positive  real  parts.  In  fact, 

VMv)  =  (2-7  t)~(*_<)BU  -  «,<;  2~fty)  -  (27«)~(*  _t)  B  (|  -  t,t;  2  7<y)  . 

For  1  <  -7  <  00,  <  is  bounded  by  0  <  t  <  [4(1  +  y/2)]  1  =  0.10355  . . .,  so  ri>1  is  well-defined. 

To  confirm  these  analytical  results,  and  to  demonstrate  of  the  accuracy  of  the  random 
choice  method,  Fig.  5  shows  a  comparison  of  the  analytical  results  with  the  results  of  solving 
the  partial  differential  equations  governing  the  rarefaction  reflection  flow.  In  this  problem, 
a  rarefaction  wave  connecting  a  pressure  po  —  1  and  density  po  =  1  to  a  vacuum  expands 
from  x  -  0  and  strikes  a  rigid  wall  at  z  =  Xo  -  1  at  time  r  =  0;  the  fluid  is  a  polytropic 
gas  with  7  =  1.4.  Two  calculations  were  performed,  one  a  free  expansion  for  which  Vo  -  0, 
and  one  in  which  the  expanding  material  has  been  shocked  to  Vo  =  ^2/7(7  -  1)cq,  which 
corresponds  to  an  infinite-strength  shock.  The  short-time  asymptotic  solution  predicts  that 
the  pressure  at  the  wall  should  rise  as  (r/tn)2'^  and  that  the  speed  of  the  reflected 
shock  should  tend  to  Ut  0.036494  •  (vn  +  2co/(7  -  1)).  In  this  computation,  the  flow  in 
the  region  0.9  <  x  <  1  and  0  -  r  <  0.3  was  simulated  using  100  mesh  cells.  Fig.  5  shows 
the  agreement  with  the  asymptotic  pressure  rise,  and  a  contour  plot  of  pressure  in  the 


Figure  5.  Pressure  at  the  wall  for  the  reflection  of  a  rarefaction  wave  expanding  into  a 
vacuum.  The  pressure  at  the  rigid  wall,  plotted  as  a  function  of  time,  is  shown  for  two 
sets  of  initial  conditions.  In  both  cases  the  initial  density  and  pressure  in  the  expanding 
material  axe  po=  1  and  p o=  1.  For  the  free  expansion,  the  initial  velocity  is  vo  =  0, 
whereas  for  the  shocked  expansion  vo  corresponds  to  driving  the  expanding  material  with 
an  infinite-strength  shock.  In  these  computations,  the  rarefaction  wave  originates  atx=0, 
and  the  flow  in  the  region  0.9  <x  <1  and  0  <t  <0.3  was  simulated  using  100  mesh 
cells.  Superimposed  are  the  short-time  asymptotic  solutions,  which  appear  as  straight 
lines.  Both  solutions  confirm  the  predicted  short-time  asymptotic  behavior  of  the  pres¬ 
sure,  viz. ,  that  it  should  rise  as  x5  when  y=  1.4.  Also,  the  solutions  are  accurate  over  8 
decades  in  pressure,  in  contrast  to  other  numerical  methods. 
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space-time  plane  shows  that  the  reflected  shock  speed  agrees  as  well.  One  may  ascertain 
that  the  short-time  asymptotic  solution  is  valid  for  r/to  <  0.1,  roughly  speaking. 

We  point  out  that  the  random  choice  solution  is  accurate  over  more  than  8  decades  in 
pressure;  this  contrasts  with  serious  discrepancies,  over  only  3  decades  in  pressure,  among 
several  numerical  methods  applied  to  this  problem,  as  reported  by  McCall.10 

Suppose  now  that  the  wall  is  replaced  with  a  thin  metallic  plate  that  is  allowed  to 
accelerate.  The  behavior  of  the  pressure  at  the  plate  is  the  same  as  for  a  rigid  wall  until 
the  plate  has  accelerated  significantly.  The  rising  pressure  generates  a  compressive  wave 
inside  the  plate,  and  we  wish  to  determine  whether  this  wave  collapses  into  a  shock  wave 
within  the  plate. 

Simple  geometry  shows  that  the  sound  characteristics  inside  the  plate  will  cross  within 
the  plate  if 

6x  =  -ii- 
de/dt 

is  less  than  the  thickness  of  the  plate.  If  the  plate  is  modeled  as  a  stiffened  gas  with 
adiabatic  exponent  70,  then 

2^0  P  +  Poo  fK  10x 

Sx= -  c  •  .  ■■■:—  (5.12) 

70-1  dp/dt 

because  the  flow  is  isentropic.  According  to  Eq.  5.3, 
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where  denotes  the  sound  speed  in  the  plate  at  zero  pressure. 

When  7  =  5/3  and  the  propellant  material  has  been  shocked  with  an  infinite-strength 
shock,  p(0)  =  51.678  po  by  the  results  above.  Therefore  if  we  consider  tungsten  (cqo  = 
0.404  cm/ps,  70  =  3.14,  pa 0  =  1  Mbar)  accelerated  by  plastic  foam  that  has  been  shocked 
to  po  =  0.24  Mbar  and  p0  =  0-12  g/cm3  and  that  is  separated  from  the  plate  by  a 
5  cm  gap,  one  calculates  that  6x  >  0.5  cm  during  the  initial  stages  of  the  acceleration, 
t  <  0.1  to  ~  0.06  ps.  This  agrees  with  the  estimate  given  by  McCall.10 


Computational  Results 

In  Figs.  6-9  are  shown  the  results  of  a  shockless  acceleration  experiment  in  which 
the  cavity  has  be  evacuated  to  10~4  bar  (the  remaining  fluid  being  a  gas  with  7  =  5/3). 
The  foam,  which  is  modeled  as  a  polytropic  gas  with  7  =  5/3,  is  initially  at  atmospheric 
pressure  and  density  0.03  g/cm3.  The  foam  is  then  shocked  to  0.24  Mbar  and  0.12  g/cm3; 
the  particle  velocity  in  the  shocked  foam  is  2.45  cm/fis. 

In  this  computation  the  tungsten  plate  is  tracked,  with  the  plate  modeled  as  a  movable 
wall  with  an  areal  mass  density  of  0.5  g/cm3,  corresponding  to  a  thickness  of  0.025  cm.  A 
3:1  mesh  refinement  is  used  for  4  <  z  <  6  cm,  with  a  total  of  300  mesh  cells:  the  mesh 
spacing  near  z  =  5  cm  is  Az  ss  0.013  cm,  while  near  z  =  10  cm  it  is  Az  «  0.035  cm. 

Figs.  6  and  7  show  the  space-time  development  of  the  flow  in  terms  of  the  contours  of 
density  and  pressure,  respectively.  At  time  t  =  0  the  shock  reaches  the  end  of  the  foam  at 
z  =  0,  whereupon  the  foam  expands  into  the  cavity.  A  very  weak  shock,  which  is  not  visible 
in  these  plots,  is  transmitted  into  the  gas  in  the  cavity  and  leads  the  expanding  foam;  its 
arrival  at  the  tungsten  plate  is  indicated  in  the  figures.  The  foam  decelerates  when  it 
strikes  the  plate,  so  a  reflected  shock  forms  and  strengthens.  The  increasing  pressure  in 
the  foam  that  builds  up  behind  the  plate  serves  to  accelerate  the  plate. 

The  pressure  behind  the  tungsten  plate  is  plotted  as  a  function  of  time  in  Fig.  8. 
Because  the  cavity  has  been  evacuated,  the  pressure  rise  is  much  smoother  than  if  no 
cavity  had  been  there.  The  velocity  of  the  tungsten  plate  is  plotted  as  a  function  of  time 
in  Fig.  9.  The  velocity  of  the  plate  at  t  =  3.85  p&  is  2.69  cm/ps. 

Fig.  8  may  be  used  to  obtain  another  estimate  for  the  crossing  distance  for  char¬ 
acteristics  in  the  plate.  The  rate  of  pressure  rise  is  estimated  from  the  graph  to  be 
dp/dt  «  1.1  Mba r/ps,  so  Eq.  5.12  gives 
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which  is  40  times  the  thickness  of  the  plate. 

In  Figs.  10-13  are  shown  the  results  of  a  shockless  acceleration  experiment  in  which 
the  cavity  is  filled  with  air  (with  7  =  5/3)  at  standard  conditions.  The  foam  is  shocked 
as  in  the  previous  computation.  In  this  computation,  both  the  foam-air  contact  and  the 
tungsten  plate  are  tracked.  The  tungsten  plate  is  modeled  as  a  movable  wall  with  an  areal 
mass  density  of  0.5  g/cm3.  A  3:1  mesh  refinement  is  used  for  4  <  z  <  6  cm,  as  before. 

Figs.  10  and  11  show  the  space-time  contours  of  density  and  pressure,  respectively. 
At  time  t  =  0  the  shock  reaches  the  end  of  the  foam  at  z  =  0,  whereupon  the  foam 
expands  into  the  cavity  and  compresses  the  air,  thereby  accelerating  the  tungsten  plate. 
Because  the  cavity  is  not  evacuated,  the  weak  shock  that  is  transmitted  to  the  air  has  more 
noticeable  effects;  its  arrival  at  the  plate  is  marked  by  a  shock  reflected  from  the  plate. 
The  contact  discontinuity  between  the  foam  and  the  air  interacts  with  waves  reflected  from 
the  plate,  which  strengthen  into  a  reflected  shock  that  overtakes  the  expansion  wave  in  the 
foam.  A  layer  of  compressed  air  moving  with  the  plate  separates  the  foam  from  the  plate. 

There  is  much  more  structure  in  the  flow  near  the  plate  in  this  computation  as  com¬ 
pared  to  the  previous  one.  The  succession  of  shock  waves  reflecting  between  the  foam-air 
contact  discontinuity  and  the  tungsten  plate  can  clearly  be  seen  by  comparing  Figs.  10 
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Figure  6.  Density  contours  for  shockJess  acceleration  using  an  evacuated  cavity.  The 
foam,  initially  at  atmospheric  pressure  and  density  0.03  g/cm3,  is  shocked  to  0.24  Mbar 
and  0.12  g/cm3;  the  particle  velocity  in  the  shocked  foam  is  2.45  cm/|is.  At  time  r  =0 
the  shock  reaches  the  end  of  the  foam  at  x  =  0,  whereupon  the  foam  expands  into  the 
cavity,  which  is  evacuated  to  10-^  bar;  the  expanding  foam  accelerates  a  tungsten  plate. 
A  reflected  shock  forms  in  the  rarefied  foam  upon  its  arrival  at  the  plate.  The  velocity  of 
the  plate  at  the  r  =  3.85  ps  is  2.69  cm/ps.  In  this  computation  the  tungsten  plate  is 
tracked,  with  the  plate  modeled  as  a  movable  wall  with  an  areal  mass  density  of  0.5 
g/cm2.  A  3:1  mesh  refinement  is  used  for  4  <x  <6  cm,  with  a  total  of  300  mesh  cells: 
the  mesh  spacing  near  x  =  5  cm  is  Ajc  =  0.013  cm,  while  near  x  -  10  cm  it  is  Ajc  =  0.035 
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Figure  7.  Pressure  contours  for  shockless  acceleration  using  an  evacuated  cavity.  The 
space-time  contours  of  pressure  are  shown  for  the  comput'tion  described  in  the  caption 
for  Fig.  6. 
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Figure  8.  Pressure  behind  the  plate.  The  pressure  behind  the  tungsten  plate  is  plotted  as 
a  function  of  time  for  the  computation  described  in  the  caption  for  Fig.  6. 
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Figure  9.  Velocity  of  the  plate.  The  velocity  of  the  tungsten  plate  is  plotted  as  a  func¬ 
tion  of  time  for  the  computation  described  in  the  caption  for  Fig.  6. 


’/  Oa 


Figure  10.  Density  contours  for  shockless  acceleration  using  an  air-filled  cavity.  The 
foam  is  shocked  as  in  the  computation  described  in  Fig.  6.  At  time  t  -  0  the  shock 
reaches  the  end  of  the  foam  at  x  =0,  whereupon  the  foam  expands  into  the  cavity  and 
compresses  the  air,  thereby  accelerating  the  tungsten  plate.  The  air  is  initially  at  standard 
conditions,  and  the  tungsten  plate  has  an  areal  mass  density  of  0.5  g/cm2.  A  weak  shock 
is  transmitted  to  the  air;  its  arrival  at  the  plate  is  marked  by  a  shock  reflected  from  the 
plate.  The  contact  discontinuity  between  the  foam  and  the  air  interacts  with  waves 
reflected  from  the  plate,  which  strengthen  into  a  reflected  shock  that  overtakes  the  expan¬ 
sion  wave  in  the  foam.  A  layer  of  compressed  air  moving  with  the  plate  separates  the 
foam  from  the  plate.  The  velocity  of  the  plate  at  the  t  =  3.85  (is  is  2.69  cm/ps.  In  this 
computation  the  foam — air  contact  and  the  tungsten  plate  are  tracked,  with  the  plate 
modeled  as  a  movable  wall.  A  3: 1  mesh  refinement  is  used  for  4  <  x  <6  cm,  as  in  Fig.  6. 
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Figure  11.  Pressure  contours  for  shockless  acceleration  using  an  air-filled  cavity.  The 
space-time  contours  of  pressure  are  shown  for  the  computation  described  in  the  caption 
for  Fig.  10.  The  succession  of  shock  waves  reflecting  between  the  foam — air  contact 
discontinuity  and  the  tungsten  plate  can  clearly  be  seen  by  comparing  this  figure  with 
Fig.  10.  Some  waves,  such  as  the  first  air  shock  that  reflects  from  the  foam — air  contact, 
are  visible  in  this  figure  but  not  in  Fig.  10,  whereas  the  secondary  contact  discontinuity  in 
the  foam  is  visible  only  in  the  density  contour  plot. 
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Figure  12.  Pressure  behind  the  plate  The  pressure  behind  the  tungsten  plate  is  plotted 
as  a  (unction  of  time  for  the  computation  described  in  the  caption  for  Fig  10  The  pres 
sure  jump  at  t  =  ()  87  is  caused  b>  the  air  shock  hitting  the  plate,  and  the  jumps  that  (<>1 
h'ss  correspond  to  the  succession  of  sh<*ck  waves  that  reflect  from  the  plate,  as  seen  m 
1  ig  1  !  The  hump  at  r~l  7  corresponds  to  the  arrival  at  the  plate  of  signals  frorn  the 
right  edge  of  the  expansion  wave  in  the  foam,  this  is  also  seen  in  Fig  1  1 
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ocity  of  the  plate  The  velocity  of  the  tungsten  plate  is  plotted  as  a  func 
the  computation  described  in  the  caption  for  Fig  10 
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a:.-!  11  Sonw  wave*,  such  u  the  first  air  shock  that  reflects  from  the  foam  air  corita<  t 
dfe  visible  in  pressure  contours  plot  but  not  in  the  density  contours  plot,  whereas  the 
s.<  "Hilary  contact  discontinuity  in  the  foam  is  visible  only  in  the  density  contour  plot 

!  tie  pressure  behind  the  tungsten  plate  is  plotted  as  a  function  of  time  in  Fig  U 
!  smooth  rise  in  pressure  seen  in  Fig.  8  has  now  become  jagged  during  the  early  stages 
'  t.<  How  development,  and  the  peak  pressure  is  higher  The  pressure  jump  at  t  0  87 
-  .  a  »>ed  by  the  air  shock  hitting  the  plate,  and  the  jumps  that  follow  correspond  to 

cession  of  shock  waves  that  reflect  from  the  plate,  as  seen  in  Fig  11  The  hump 

at  f  1  7  corresponds  to  the  arrival  at  the  plate  of  signals  from  the  right  edge  of  the 
expansion  wave  ir  the  foam:  this  is  also  seen  in  Pig  11  Pig  13  shows  the  velocity  of  the 
tungsten  plate,  plotted  as  a  function  of  time  The  kinks  in  this  plot  correspond  to  the 
Tm  ontmuities  in  P  ig  12  The  velocity  of  the  plate  at  t  3  85  jcs  is  again  2.69  cm  ^s 

Although  the  short  time  development  of  the  flow  changes  markedly  when  the  cavity 
pressure  is  raised  to  as  small  a  pressure  as  10  6  Mbar,  the  long-time  behavior  of  the  flow 
and  the  accelerating  plate  is  little  affected  This  is  evident  in  the  contour  plots  and  the 
pmts  of  the  velocity  of  and  pressure  behind  the  plate.  In  fact,  at  late  times  the  velocity 
and  pressure  for  the  air-filled  cavity  agree  quantitatively  with  these  quantities  when  the 
.  .iv  'v  is  evacuated 

!  He  reverberating  shock  waves  that  are  present  in  the  flow  when  the  cavity  has  not 

•  -  -  i;  i*v a<  nated  may  cause  damage  to  the  accelerating  plate,  for  instance  through  shock 

'  -a'.ng  I'o  -v  certain  the  effect  that  these  air  shocks  have  on  the  plate,  the  experiment 
W  tr.e  air  tided  cavity  was  simulated,  during  the  early  stages  of  the  interaction,  with 

•  ’  •  •  mgsten  plate  modeled  as  a  thin  (0.025  cm)  layer  of  fluid  In  this  computation  the 

.at  mi,  of  stat«‘  of  the  tungsten  is  a  stiffened  gas,  with  3  14  and  p  x  1.0  Mbar.  The 
Tw"  sides  of  the  plate  are  tracked,  as  is  the  foam-air  contact  discontinuity.  A  10:1  mesh 
o  '.iiement  is  used  for  4  95  *  x  <  5  25  cm,  with  a  total  of  300  mesh  cells.  With  this  mesh 
r<  K.ement  there  are  approximately  10  mesh  cells  inside  of  the  plate. 

1  he  results  are  shown  in  Figs  14-18.  Fig.  14  shows  the  space-time  contours  of  pressure, 

•  which  the  shocks  that  form  inside  the  plate  are  evident.  An  expanded  view  of  the 
pre^ure  contours  in  the  vicinity  of  the  accelerating  plate  is  shown  in  Fig.  15.  The  detailed 
.:.*»•  r r. a!  structure  of  the  How  inside  the  plate  is  evident  in  this  plot.  The  corresponding 
; . . .  it  of  the  pressure  behind  the  plate,  plotted  as  a  function  of  time,  is  shown  in  Fig.  16. 

!  pressure  jump  at  t  0  87  is  caused  by  the  air  shock  hitting  the  plate,  and  the  jumps 
Wat  foiiow  correspond  to  the  succession  of  shock  waves  that  reflect  from  the  plate,  as  seen 
i  g^-  H  and  1 5 

1  n*'  relative  velocity  of  the  two  sides  of  the  plate  is  plotted  as  a  function  of  time  in 
Kg  i  7  \  positive  relative  velocity  indicates  expansion  of  the  plate  The  jumps  in  the 
i  >  pariMor.  velocity  correspond  to  waves  in  the  air  arriving  at  the  plate  and  to  waves  in 

•  e  ;i,a'*’  reflecting  from  the  ends  of  the  plate,  as  seen  in  Fig.  15.  Shock  waves  in  the  air 
‘na’  v w i k ••  the  bai  k  end  of  the  plate  transmit  a  shock  wave  into  the  plate  and  cause  the 
•.  a ' :  v  veioc  it y  to  drop  The  transmitted  shock  reflects  from  the  front  end  of  the  plate  as 
a  ran  ’ ac  lion  wave  and  causes  the  relative  velocity  to  rise.  This  simple  ringing  of  the  plate 

’  r,  ,.'s  ■  i r, 1 1 1  mtera<  turns  between  waves  creates  broad  rarefaction  waves  that  complicate 
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Figure  14.  Pressure  contours  for  shockless  acceleration  using  an  air-filled  cavity.  The 
space-time  contours  of  pressure  are  shown  for  the  early-iime  development  in  the  compu¬ 
tation  described  in  the  caption  for  Fig.  10.  In  this  computation  the  tungsten  plate  is 
modeled  as  a  stiffened  gas.  The  two  sides  of  the  plate  are  tracked,  as  is  the  foam — air 
contact  discontinuity.  A  10:1  mesh  refinement  is  used  for  4.95  <  x  <5.45  cm,  with  a 
total  of  300  mesh  cells.  With  this  mesh  refinement  there  are  approximately  10  mesh  cells 
inside  of  the  plate.  The  combination  of  mesh  refinement,  front  tracking,  and  the  random 
choice  method  is  crucial  to  resolving  the  shock  waves  that  can  be  seen  propagating  inside 
the  plate. 
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Figure  15.  Expanded  view  of  the  plate.  The  pressure  contours  for  the  computjtioi 
described  in  the  caption  for  Fig.  14  are  shown  in  the  vicinity  of  the  accelerating  plate 
The  detailed  internal  structure  of  the  flow  inside  the  plate  is  evident  in  this  plot 
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Figure  16.  Pressure  behind  the  plate.  The  pressure  behind  the  tungsten  plate  is  plotted 
as  a  function  of  time  for  the  computation  described  in  the  caption  for  Fig.  14.  The  pres- 
suR-  jump  at  t  =  0  87  is  caused  by  the  air  shock  hitting  the  plate,  and  the  jumps  that  fol¬ 
low  correspond  to  the  succession  of  shock  waves  that  reflect  from  the  plate,  as  seen  in 
1  ics  14  and  15 


Figure  17.  Expansion  of  the  plate.  The  relative  velocity  of  the  two  sides  of  the  plate  is 
plotted  as  a  function  of  time  for  the  computation  described  in  the  caption  for  Fig.  14.  A 
positive  relative  velocity  indicates  expansion  of  the  plate.  The  jumps  in  the  expansion 
velocity  correspond  to  waves  in  the  air  arriving  at  the  plate  and  to  waves  in  the  plate 
reflecting  from  the  ends  of  the  plate,  as  seen  in  Fig.  15. 
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Figure  18.  Shock  heating  inside  the  tungsten  plate.  An  approximation  A Q  to  the  heat 
deposition  inside  the  plate  is  shown.  The  maximum  and  the  spatial  average  of  A Q  across 
the  plate  is  plotted  as  a  function  of  time.  For  comparison,  the  latent  heat  of  fusion  of 
tungsten  at  standard  conditions  is  0.0019  cm2/us2. 


Although  the  ringing  of  the  plate  is  rather  violent,  the  rarefaction  waves  never  combine 
to  put  the  material  in  tension,  as  they  do  in  a  spall  experiment.  (This  is  confirmed  by 
plotting  the  minimum  pressure  in  the  plate  as  a  function  of  time.)  Within  the  confines  of 
this  simplistic  model  for  the  equation  of  state,  therefore,  it  is  expected  that  the  plate  will 
not  spall  when  it  is  accelerated  in  the  manner  being  investigated. 

As  another  diagnostic  measure  for  possible  structural  damage  to  the  plate,  the  heat 
deposition  in  the  plate  due  to  shock  waves  should  be  determined.  This  heat  deposition  is  to 
be  compared  to  the  heat  of  melting  of  tungsten  to  determine  whether  the  plate  melts.  Since 
the  plate  is  being  accelerated  by  a  light  fluid,  it  is  subject  to  Rayleigh-Taylor  instability. 
Were  the  plate  to  melt,  its  material  strength  might  not  sufficiently  dampen  the  growth  of 
these  instabilities,  and  it  could  be  destroyed. 

The  stiffened  gas  equation  of  state  is  not  complete,  in  that  the  temperature  and 
entropy  is  not  determined.  In  place  of  using  a  complete  equation  of  state,  an  approximation 
to  the  heat  deposition  in  the  plate  is  determined  by  comparing  the  computed  specific 
energy  in  the  plate  with  the  specific  energy  that  would  be  in  the  plate  if  the  plate  had 
been  isentropically  compressed.  Let 
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7  ~  1  P 

be  the  specific  energy  in  a  given  computational  mesh  cell  inside  the  plate,  and  let 

1  (p/PO^PO  +  Poo)  +  (7  -  l)Poo 


Es0  = 


7  -  1 


be  the  corresponding  specific  energy  of  the  cell  were  it  to  have  been  compressed  isen¬ 
tropically  from  its  initial  state  with  density  po  and  pressure  po-  (See  Eq.  2.1.)  We  then 
define 

AQ  =  E-  ESo  «  TAS 

to  be  an  effective  heat  deposition  in  the  cell.  The  maximum  of  this  quantity  over  the  plate 
and  its  spatial  average  is  plotted  as  functions  of  time  in  Fig.  18.  Since  the  latent  heat  of 
fusion  of  tungsten  at  standard  conditions  is  0.0019  cm2//is2,  which  is  of  the  same  order  of 
magnitude  as  AQ,  no  definitive  conclusion  can  be  drawn  from  the  calculation,  although  it 
suggests  that  the  plate  might  melt. 


Conclusions 

These  computations  demonstrate  the  ability  of  our  numerical  method  to  resolve  the 
highly  complicated  flows  that  arise  in  McCall’s  method  for  accelerating  thin  plates.  The 
combination  of  mesh  refinement,  front  tracking,  and  the  random  choice  method  is  crucial 
to  resolving  shock  waves  propagating  inside  the  plate. 

The  short-time  asymptotic  solution  and  the  computational  results  confirm  that  shock 
waves  do  not  form  inside  a  sufficiently  thin  plate,  provided  that  the  cavity  between  the 
propellant  foam  and  the  plate  has  been  evacuated.  The  computations  do  show,  however, 
that  it  is  critical  to  have  a  good  vacuum  in  the  cavity:  pressures  of  only  10-6  on  the  scale 
of  pressures  being  considered  give  rise  to  significant  shock  waves  in  the  plate.  These  shock 
waves  cause  violent  ringing  of  the  plate,  but  they  never  put  the  plate  into  tension;  thus 
according  to  the  model  used,  the  plate  does  not  spall.  The  shock  waves  also  heat  the  plate, 
and  a  simple  estimate  of  this  effect  indicates  that  the  plate  might  melt.  This  melting  would 
make  the  plate  more  susceptible  to  Rayleigh-Taylor  instability. 

In  combination  with  a  better  model  of  the  structural  properties  of  the  tungsten  plate, 
the  results  of  these  hydrodynamic  calculations  can  be  used  to  more  conclusively  determine 
the  feasibility  of  McCall’s  method  for  shockless  acceleration  of  thin  plates. 
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